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EVOLUTION OF THE CLUSTER X-RAY SCALING RELATIONS SINCE Z > 0.4 
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ABSTRACT 

We derive correlations between X-ray temperature, luminosity, and gas mass for a sample of 22 distant, z > 
0.4, galaxy clusters observed with Chandra. We detect evolution in all three correlations between z > 0.4 and 
the present epoch. In particular, in the D. = 0.3, A = 0.7 cosmology, the luminosity corresponding to a fixed 
temperature scales approximately as (1 +z) 15±0 - 3 ; the gas mass for a fixed luminosity scales as (1 +z) , 
and the gas mass for a fixed temperature scales as (1 +z) — °' 5±0 ' 4 (all uncertainties are 90% confidence). We briefly 
discuss the implication of these results for cluster evolution models. 

Subject headings: galaxies: clusters: general — surveys — X-rays: galaxies 



1. INTRODUCTION 

Correlations between the X-ray properties of galaxy clus- 
ters are useful statistical tools which allow to study the clus- 
ter physics. The best-studied correlation for the low-redshift 
clusters is that of the X-ray luminosity with the temperature 
(e.g., Mushotzky 1984, David et al. 1993, Markevitch 1998, 
Arnaud & Evrard 1999). The mass of the intracluster gas corre- 
lates with both temperature (Mohr, Mathiesen & Evrard 1999, 
Vikhlinin, Forman & Jones 1999) and X-ray luminosity (Vo- 
evodkin, Vikhlinin & Pavlinsky 2002b). The observed tightness 
of these correlations indicates a similar formation history for all 
clusters and is consistent with the predictions of the self-similar 
models of the cluster formation. 

However, the details of these scaling relations are different 
from the predictions of the self-similar theory. The best-known 
example is the slope of the L — T correlation: it is observed that 
L °< T 2n for hot clusters (e.g., Markevitch 1998), while the- 
ory predicts L °< T 2 (e.g., Kaiser 1981). Such deviations may 
be due to non-gravitational processes such as preheating (e.g., 
Cavaliere, Menci & Tozzi 1997) or radiative cooling and feed- 
back from star formation (Voit & Bryan 2001). Observing the 
evolution of scaling relations can provide useful constraints on 
such models. The scaling relations at high redshift also are of 
great value for cosmological studies based on cluster evolution. 
They provide the means to convert an easily observed X-ray 
luminosity function into the more cosmologically useful tem- 
perature or mass functions (e.g., Borgani et al. 2001). 

Most of the previous high redshift studies have focused on 
the L — T relation. Mushotzky & Scharf (1997) analyzed a large 
sample of distant clusters observed with ASCA (most at z ~ 0.3, 
with a few at z > 0.4) and found no evidence for evolution in 
the L—T relation. Several results for a small number of dis- 
tant cluster observed by Chandra have been published recently. 
Borgani et al. (2001) analyzed a sample of 7 clusters at z > 0.5 
and concluded that the data allow at most a very mild evolution 
— if the luminosity for the given temperature is L(z) « ( 1 +z) A , 
then A < 1 . A similar conclusion has been reached by Holden 
et al. (2002) from an analysis of 12 clusters at z > 0.7. 

The L — T relation at both low and high redshifts has a large 
intrinsic scatter, comparable to the expected evolutionary ef- 
fects. The scatter in the low-redshift L — T relation is signifi- 



cantly reduced when the central cooling regions of the clusters 
are excised from both the luminosity and temperature measure- 
ments (Fabian et al. 1994, Markevitch 1998). Therefore, it is 
desirable to exclude the cooling cores in the distant clusters 
because this too may reduce the scatter and thus more easily 
expose any evolution. This task is feasible only with the Chan- 
dra's arcsecond angular resolution. 

As of Spring 2002, Chandra had observed 22 clusters at 
z > 0.4 with sufficient exposure for accurate temperature mea- 
surements. Most of this sample is derived from flux-limited 
X-ray surveys: 7 clusters from the EMSS (Henry et al. 1992), 
1 1 from the ROSAT serendipitous surveys, 160 deg 2 (Vikhlinin 
et al. 1998), RDCS (Rosati et al. 1998), WARPS (Ebeling et al. 
2000), and 2 from the ROSAT All-Sky Survey. We use these 
Chandra observations to derive correlations between the X-ray 
luminosity, temperature, and gas mass at z > 0.4. 

We use Hq = 50 km s _1 Mpc~' throughout. 

2. DATA ANALYSIS 

Chandra data were reduced in a standard manner (see Marke- 
vitch & Vikhlinin 2001 for a fuller description of our proce- 
dures). Spectral analysis was performed in the 0.8-10 keV 
band, to minimize the calibration uncertainties due to time- 
dependent decrease in the low energy quantum efficiency 3 . The 
effective area below 1.8 keV in the front-illuminated CCDs was 
corrected by an empirical factor of 0.93 to improve the cross- 
calibration with the back-illuminated CCDs (see Markevitch & 
Vikhlinin 2001 for details). Background maps were generated 
using the prescriptions of Markevitch (2002). Cluster spec- 
tra were fit with the MEKAL model, with absorption fixed at 
the Galactic value, and abundance fixed at 0.3 Solar unless the 
quality of the data allowed a direct measurement. 

The imaging analysis was performed in the 0.7-2 keV energy 
band, but we also checked 2-7 ke V images to excise hard, self- 
absorbed sources which may contaminate the spectral analysis 
of the faintest clusters. 

Our aim was to compare the high-redshift L —T relation with 
the low-redshift one from Markevitch (1998), and so we fol- 
lowed this analysis as closely as possible. Temperatures were 
measured by fitting a spectrum integrated within a radius of 
0.5-1 Mpc with a single-temperature MEKAL model; this re- 
gion contained at least 75% of the total cluster flux. The bolo- 
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Table 1 
Cluster sample 



Name 


z 


T 


^0 5 2 


L bol 


fa 






(keV) 






(lO r4 M ) 


MS 0016+1609 


0.541 


9.9±0.5 


22.8 


113.3 


6 43 ±0 65 


MS 0302+1658 


0.424 


3.6 ± 0.5 


4.7 


10.6 


i 07 + 40 


MS 0451-0305 


0.537 


8.1 ±0.8 


20.7 


91.7 


3 fj« + n 77 


MS 1054—0321 


0.823 


7.8±0.6 


16.5 


70.9 


9 5R + 37 


MS 1 1 37+6625 


0.782 


6.3 ±0.4 


8.4 


32.4 


1 41 +0 28 


MS 1621+2640 


0.426 


7.6±0.9 


6.3 


27.0 


2 89 + 62 


MS 2053-0449 


0.583 


5.2 ±0.7 


3.5 


10.8 


n q<; + 32 


Tl 1 1 20+2326 


0.562 


4.8±0.5 


3.7 


12.5 


1.19±0.27 


rj 1221 +491 X 

VI.. 1 1 T*Ty 1 □ 


0.700 


7.2±0.6 


7.0 


28.7 


2 oi +n 36 


CL 1416+4446 


0.400 


3.7±0.3 


4.2 


8.9 


1.42 ±0.38 


TT 1 524+0957 


0.516 


5.1 ±0.6 


4.5 


15.7 


1 fi7 + n 40 


CL 1701+6421 


0.453 


5.8±0.5 


4.9 


15.9 


1.81 ±0.47 


CL 0848+4456 


0.574 


2.7±0.3 


10.6 


38.8 


0.36±0.13 


WARPS 0152-1357 


0.833 


5.8±0.6 


1.2 


3.1 


2.88±0.55 


RDCS 0848+4452 


1.261 


4.7 ±1.0 


1.8 


6.0 


0.20 ±0.08 


RDCS 0910+5422 


1.100 


3.5 ±0.7 


2.0 


5.9 


0.26 ±0.11 


RDCS 1317+2911 


0.805 


2.2±0.5 


0.8 


2.0 


0.21 ±0.09 


RDCS 1350.0+6007 


0.805 


4.3±0.6 


4.2 


13.2 


1.04 ±0.33 


RASS 1347-114 


0.451 


14.1 ±0.9 


60.1 


260.4 


8.77 ±1.60 


RASS 1716+6708 


0.813 


6.6±0.8 


7.2 


28.8 


1.25 ±0.33 


3C295 


0.460 


5.3 ±0.5 


9.1 


16.3 


1.51 ±0.48 


CL0024+17 


0.394 


4.8±0.6 


3.1 


9.2 


1.24 ±0.37 



° — Total X-ray luminosity in the 0.5-2 keV band within the 2 Mpc radius, 10 44 erg s * — Bolometric luminosity within the 
2 Mpc radius excluding the central cooling regions, 10 44 erg s _1 . All quantities are computed for the Q. — 0.3, A = 0.7 cosmology. 



metric X-ray luminosity was computed by normalizing the best- 
fit spectral model by the observed total 0.7-2 keV count rate 
within 2 Mpc. The extrapolation of a typical j3 -model profile 
shows that at most 4% of the X-ray flux originates outside this 
radius, which we ignore. The clusters with sharply peaked sur- 
face brightness profiles are likely to contain cooling cores, and 
so we excluded the central 100 kpc regions in such clusters from 
both the spectral analysis and count rate computations, multi- 
plying the measured count rates by 1 .06 to account for the flux 
within r < 100 kpc in a typical j3 -model cluster. 




Chandra T, keV 

FIG. 1. — Comparison of Chandra and ASCA temperatures. ASCA cool- 
ing flow corrected values for A401, A478, A644, A754, A85, A780, A2029, 
A2256, A2597, A3376, A3558, A3667, and MKW3s are from Markevitch 
(1998; circles). Triangles represent non- or weak cooling flow clusters A1060, 
A2147, A2218, A2255, and AWM7 with the ASCA values from White (2000) 
without any cooling flow correction. Uncertainties are 90% confidence. 



It is important to verify that there is no systematic difference 
between our Chandra and earlier ASCA temperature measure- 
ments. Using the procedure outlined above, we derived temper- 
atures for a number of nearby clusters from the ASCA samples 
of Markevitch (1998) and White (2000) and confirmed that on 
average, the agreement between Chandra and ASCA measure- 
ments is at a level of 5% or better (Fig. 1). 

3. RESULTS 

The resulting L — T correlation in the £2 = 0.3, A = 0.7 cos- 
mology is shown in Fig. 2. The scatter in the correlation is 
very small — it is almost consistent with the measurement un- 
certainties. The slope is consistent with the low-redshift re- 
lation, but the normalization is significantly different — for 
the same temperature, clusters at high redshift are more lumi- 
nous. If the z-dependence of the correlation is parameterized 
as L °c (1 +z) Alt T a with a fixed at its local value of 2.64, 
then A LT = 1.5 ±0.3 (90% confidence) for this cosmology. For 
Q = 1, A = 0, evolution of the L - T relation is weaker, but still 
significant (Table 2). 

Let us now consider correlations involving the gas mass. As- 
suming spherical symmetry, the radial profile of the gas mass 
in a cluster is trivially derived from the observed X-ray sur- 
face brightness profile. The main question is: within which 
radius to measure the gas mass? We chose the radius defined 
by the mean gas overdensity over the average baryon density 
in the Universe 8 S = M g (r)/ (4-/3n r 3 (p h )) = 324 (/z/0.5) 1 / 2 , 
where M g (r) is the gas mass within the radius r, and (p b ) = 
5.55M Q kpc~ 3 x (1 +z) 3 is the average baryon density from 
Big Bang nucleosynthesis theory (Buries, Nollett, & Turner 
2001). This choice is motivated by two considerations. First, if 
the mass fraction of baryons in clusters is representative of that 
for the Universe as a whole, then at large distances from the 
cluster center, 8 g should equal the mean total mass overdensity, 
8 m , and at the same time 8 g is easily determined from the X- 
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FIG. 2. — Correlation of the cluster temperature and bolometric luminosity. Both quantities are measured excluding any central cooling regions of the clusters. 
Filled circles correspond to clusters with 0.4 < z < 0.7, and open circles to clusters at z > 0.7. The solid line shows the L — T correlation for the low-redshift clusters 
(Markevitch 1998). The right panel shows the effect of removing the best-fit evolution given in Table 2. 



ray imaging data (our procedure is described in Vikhlinin et al. 
1999). Strictly speaking, the baryon mass should include the 
stellar material but we ignore this minor contributor. Second, 
Jenkins et al. (2001) developed a universal model for the clus- 
ter total mass function for masses corresponding to 8,„ = 324 at 
the redshift of observation, therefore our values can be readily 
compared to theoretical predictions. 




T, keV 

FIG. 3. — Correlation of cluster temperature and gas mass within the radius 
of mean overdensity of 324. Symbols same as in Fig. 2. Solid line shows the 
M gas - T correlation for the low-redshift clusters (Voevodkin et al. 2002a). 

Figure 3 shows the correlation of M g324 with the cluster 
temperature. The solid line shows this relation for a sam- 
ple of low-redshift clusters (Voevodkin et al. 2002a). For the 
D. = 0.3, A = 0.7 cosmology, the evolution is very weak. If 
the z-dependence in the correlation is parameterized as M g °c 
(1 +z) a mttI } with j3 fixed on its local value of 1.76, then 



A MT = 0.5 ±0.4. For Q = 1, A = 0, the derived gas masses 
are significantly smaller and there is a strong evolution in the 
M g - T relation, A MT = 1 .5 ± 0.4. 

Figure 4 shows a correlation of M 324 with the X-ray lu- 
minosity. This relation is potentially useful for estimating the 
mass function from the X-ray luminosity function for distant 
clusters. Therefore, we used in this case the total cluster lu- 
minosity (i.e. not excluding the central regions) in the 0.5- 
2 keV energy band — the quantity commonly measured in the 
serendipitous ROSAT cluster surveys. The solid line in Fig. 4 
shows the M — L tot correlation for low-redshift clusters (Vo- 
evodkin et al. 2002b). The M g — L relation shows a strong evo- 
lution in all cosmologies. If its z-dependence is parameterized 
as a power law, M °< (1+ z) A " L L r with y fixed at its local value 
of 0.83, then A ML w 2 ±0.4 (Table 2). 

Table 2 

Evolution of scaling relations in different 
cosmologies 



Cosmology 


% 


A LT 


A MT 


A ML 


£2 = 0.3, A = 0.7 


-0.55 


1.5 ±0.3 


0.5 ±0.4 


1.8±0.4 


a = 1,A = 


0.5 


0.6 ±0.3 


1.5 ±0.4 


2.1 ±0.4 



For other cosmologies, the power law slopes in the relations L oc ( 1 + z) Alt T a , 
, M g oc ( 1 + z) ml V can be approximated by interpolating 
over q . Uncertainties are at the 90% confidence level for a single parameter. 



4. DISCUSSION AND CONCLUSIONS 

We have used Chandra observations of a sample of 22 dis- 
tant, z > 0.4, clusters to show that the correlations between the 
cluster temperature, luminosity, and gas mass evolve signifi- 
cantly with respect to the low-redshift relations. 

Our detection of significant evolution in the L—T relation 
appears to contradict some other recent Chandra studies (Bor- 
gani et al. 2001, Holden et al. 2002). The difference between 
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10" io 45 

L x [0.5-2 keV], erg s" 




10" 10 45 

L x [0.5-2 keV], erg s" 



FlG. 4. — Correlation of cluster luminosity in the 0.5-2 keV band and gas mass within the radius of mean overdensity of 324. Symbols same as in Fig. 2. The 
luminosities include the central cooling regions. Solid line shows the M — L correlation for low-redshift clusters, and dotted lines correspond to a 1-0" scatter about 
the mean relation (Voevodkin et al. 2002b). The right panel shows the effect of removing the best-fit evolution given in Table 2. 



their and our results can be traced mostly to a more consis- 
tent comparison of the high- and low-redshift samples, such as 
exclusion of the cool cores and extraction of the luminosities 
in the 2 Mpc aperture. Note that the results of Novicki et al. 
(2002) who self-consistently use ASCA data for the nearby and 
distant clusters, agree well with our L — T evolution. 

It is theoretically expected that within clusters, the baryon 
contribution to the total mass, f b , should be close to the aver- 
age value in the Universe (e.g., White et al. 1993). This notion 
continues to gain observational support (most recently, Allen, 
Schmidt & Fabian 2002). If f h is indeed constant, the evolu- 
tions in the M — L and M — T relations involving the gas mass 
or total mass should be identical. 

The observed evolution of the cluster M g — T — L correla- 
tions indicates that clusters at high redshift were systematically 
denser than at present — hotter and more luminous for a given 
mass, as expected in a theory of the hierarchical self-similar for- 
mation. However, the details of the observed evolution contra- 
dict the self-similar predictions. For example, the standard the- 



ory (e.g., Bryan & Norman 1998) predicts that for a given tem- 
perature, the product H(z)M A (T,z) should be constant, where 
M A is the mass measured within the radius of the overdensity A 
with respect to the critical density at redshift z, and H(z) is the 
Hubble constant. For the realistic cluster density profiles, this 
implies that approximately M g (z, T)<x(l +z)~ 3 / 2 in almost any 
cosmology, where Mg corresponds to the mean overdensity 8 
relative to the average density at redshift z- However, we ob- 
serve little evolution in the M g 5 — T relation for the currently 
favored £2 = 0.3, A = 0.7 cosmology, seemingly at odds with 
the theoretical prediction. This possibly indicates the impor- 
tance of non-gravitational processes for heating the intracluster 
gas that would change T but are unlikely to modify M g , just as 
is inferred from the scaling relations for low-redshift clusters. 

This work was supported by NASA grant NAG5-9217 and 
contract NAS8-39073. We thank J. P. Henry and M. Novicki 
for sharing their results prior to publication. 
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